MPS-RR 2001-24 



A String Bit Hamiltonian Approach 

to 

Two-Dimensional Quantum Gravity 

B. Durhuus'^'^fl and C.-W. H. Lee^'^fl 

" MaPhySto — Centre of Mathematical Physics and Stochastic^ 
^ Department of Mathematics, University of Copenhagen, 
Universitetsparken 5, DK-2100 Copenhagen, Denmark. 

August 21, 2001 
Abstract 

Motivated by the formalism of string bit models, or quantum matrix models, 
we study a class of simple Hamiltonian models of quantum gravity type in two 
space-time dimensions. These string bit models are special cases of a more 
abstract class of models defined in terms of the SI2 subalgebra of the Virasoro 
algebra. They turn out to be solvable and their scaling limit coincides in 
special cases with known transfer matrix models of two-dimensional quantum 
gravity. 



PACS numbers: 04.60.Kz, 04.60.Nc. 

Keywords: quantum gravity, string bit model, large- limit, scaling limit, 
s/2 algebra. 



^e-mail address: durhuus@math.ku.dk 
^e-mail address: lee@math.ku.dk 

■'funded by a grant from the Danish National Research Foundation. 



1 



1 Introduction 



So far the most productive approach to two-dimensional quantum gravity 
has been in terms of path integrals. Specifically, the (Euclidean) path inte- 
gral formulation of Liouville gravity together with its regularized versions 
involving random triangulations, or the equivalent matrix models (see, e.g., 
Ref. 0), have permitted detailed analysis of two-dimensional quantum grav- 
ity coupled to matter fields with central charge c < 1, in particular the loop 
correlation functions (or Hartle-Hawking wave functions) and, in some cases, 
the fractal characterstics of space-time. Importantly, and as a rare instance 
in quantum field theory, the discrete models have provided analytic insights 
that presently seem out of reach of continuum methods, the most striking 
of which is perhaps the double scaling limit as a procedure to incorporate 
space-times of arbitrary topology. (See Ref.0] and the references therein.) 

One of the main motivations for studying quantum gravity in two dimen- 
sions is based on the hope that it may serve as a testing ground for ideas and 
methods extendable to higher dimensions. It is, indeed, straightforward to 
set up discrete models of quantum gravity in arbitrary dimensions in terms 
of random triangulations, but up to now very few analytic results have been 
obtained and even very basic questions are left open. A number of numeri- 
cal investigations have, however, been carried out. See, e.g., Ref.0] and the 
references therein. 

The question arises, naturally, if there exist alternative formulations of 
the regularized models, or closely related models, that are better tailored for 
generalizations. A second, and independent, purpose of such reformulations 
would be to make comparisons possible with continuum approaches other 
than the path integral quantization, in particular canonical quantization Q . 
In this paper we address this question in two dimensions and introduce a class 
of (Euclidean) Hamiltonian models of regularized two-dimensional quantum 
gravity. We do not claim to resolve the above mentioned questions, but we 
will show that the proposed models provide a Hamiltonian alternative to 
the discrete path integral (or transfer matrix) approach to a class of mod- 
els introduced recently under the name Lorentzian gravity 0, |^, which on 
the other hand are closely related to the full randomly triangulated models 
mentioned at the beginning. More precisely, we will show that the contin- 
uum limits of the Lorentzian models can be obtained from our Hamiltonian 
models. Indeed, we will define and solve a more abstract class of models and 
introduce a formalism opening up the way for even further generalizations. 
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Detailed analysis of these generalized models as well as generalizations to 
higher dimensions is, however, outside the scope of the present paper. 

Here is the basic idea of this work. Recall that the primary goal of the 
Hamiltonian formulation of quantum gravity is to account for the "time"- 
evolution of a space-like universe of fixed topology^. Restricting first to con- 
nected and compact space-like universes implies that the two-dimensional 
space-time has either the topology of a strip, corresponding to equal-time 
slices that are open line segments, or a cylinder with circles as equal-time 
slices. Since the only reparametrization-invariant quantity defined by a met- 
ric on a one- dimensional (connected) manifold is its volume, a natural way 
to discretise the spacial metric degree of freedom is to introduce a distance 
cutoff a > and consider the equal-time slices to be polygonal lines or loops, 
respectively, with volume n ■ a, where the integer n is the number of links 
in the slice. Keeping a fixed, we associate with each such space-like universe 
of volume n ■ a a, pure quantum state | n) , and these are assumed to form an 
orthogonal basis for the Hilbert space T of states. The Hamiltonian acting 
on T is chosen in such a way that the action couples adjacent links only. 

It turns out that models of this type may conveniently be generated by 
a special variant of the string-bit formalism, whose basic ingredients are 
annihilation and creation operators that can be interpreted as annihilating 
or creating links in the equal-time slices. (We will explain the necessary 
details in Sections || and 0.) String bit models were originally developed as a 
means of regularizing string theory [7|. They provide suitable frameworks 
for quantum chromodynamics and quantum spin chain models 0, too. 
Our variant may be envisaged as generalized quantum spin chain models in 
which the numbers of "spins", i.e. links, are variable. In this sense, the 
relationship between Lorentzian gravity models and those string bit models 
equivalent to them is analogous to that between the six-vertex and the XXZ 
model ; a better understanding of one class of models will spur the study 



of the other. 

The paper is organised as follows. In Section 2, the simplest possible 
Hamiltonian model in the case when space-time is a strip will be solved in 
the continuum limit. This model turns out to coincide with the correspond- 
ing Lorentzian gravity model considered in 0; in Section ^, we will consider 
a Hamiltonian model for cylindrical space-time which is not spacially homo- 
geneous (or cyclically symmetric). The model will be shown to reproduce, in 

"^Note that the models considered in this paper are all within the Euclidean framework. 
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the continuum limit, the Lorentzian model with a marked link in the initial 
space-like slice considered in in Section 4, we will consider the cyclically 
symmetric version of the previously mentioned model and show that this, as 
well as the model in Section 2, can be obtained as special cases of a more gen- 
eral class of models expressed in terms of the s/2-generators of the Virasoro 
algebra in a certain class of highest weight representations, the Hamiltonian 
being of the form 

H = Lo + \Li + AL_i. (1) 

In Section 5, we will solve this model. In particular, both the two-loop am- 
plitude of the continuum Liouville gravity model in Ref. [|ll| and the so-called 
p-seamed correlation functions of Ref. will be obtained as special cases; fi- 
nally, we will discuss briefly further generalisations and future developments 
in Section |^. 



2 Space-time with Boundaries 

In this section, we are going to consider the quantum mechanics of a space- 
time with the topology of a strip whose Hamiltonian is given by 



Hi 
2 



Tr 



a' a + 



N 



N 



(2) 



for = 00. Before explicating the various terms in Eq.(^), let us briefly 



review what string bit models are. We will largely follow Refs.||12|| and |T3 
•with a few modifications in definitions and notations. 

Consider an x matrix of creation operators. Its matrix entry is 
written as aj^'^^, where /^i and /i2 are row and column indices, respectively, 
and can take any integer values between 1 and inclusive. The Hermitian 
conjugate of this matrix is an A^ x A^ matrix of annihilation operators whose 
matrix entries are written in the form a^^. The creation and annihilation 
operators satisfy the canonical commutation relations 



^2 ' A'4 



/i4 /i2 



All other commutators involving these operators vanish. 

In addition, consider a 1 x A^ row vector and an A^ x 1 column vector of 
creation operators. Their components are written as q^^ and respectively. 
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Their Hermitian conjugates are an x 1 column vector and a 1 x row 
vector of annihilation operators. Their components take the form and q^, 
respectively. These operators satisfy the canonical commutation relations 



and 



"mi 



^^2 



All other commutators involving them, including those involving a^^^ or a^^ 
as well, vanish. 

Let \VL) be a vacuum state annihilated by all annihilation operators, and 
define Ti as the Hilbert space spanned by all states of the form 



\n) 



1 



jV(n+l)/2 
1 



Ml M2 



Mn I ' 



gt(at)Vlf^), 



Ar('^+i)/2 

where n is a positive integer and the summation convention is adopted for 
all row and column indices. The choice of notation I ) i for vectors in Ti will 

2 2 

be explained in Section 4. The inner product (■ | ■)! on Ti is fixed uniquely 
by the commutation relations and = 1. Since |T2 



lim (m|n) 



the states |1) i, |2) i, . . . , and so on form an orthonormal basis of Ti in the 

2 2 2 

large- A^ limit, which is the limit we are considering. We think of |n) i as the 
quantum state of a universe made up of n interior links and two boundary 
links. 

There are various kinds of natural operators acting on Ti. We define 



1 



]\f{k+e-2)/2 Ml M2 

I Tr \(a^)''a'^ 



Mfc I'l ^l-i 



7M1 



(3) 



where k and ^ are any positive integers, and Tr denotes the trace in index 
space (and not on Ti). Moreover, we let 



and 



'0 



q^q 



(4) 
(5) 
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where the superscript t denotes transposition. In the large- hmit T2\, 









fo 




n) 1 








2 






'■0 


n) 1 

' 2 




n) 1 

' 2 


oo 


n) 1 







if £ > n; 
if £ < n; 



(6) 

where (9(1/A^) consists of terms whose norms and whose inner products with 
any |n)i are of the order of at most These terms are thus neghgible in 

the large- limit. We can see from Eq.(|^) that a\ replaces any I adjacent in- 
terior links with k adjacent interior links and annihilates |n) i if £ > n. /q and 
rg annihilate the left and right boundary links, respectively, and then create 
them back. Hence, both /{] and rfj effectively act as the identity operator on 
T^, but we will nevertheless display them explicitly below. 

Using Eqs.d), (|), and d), we can rewrite the Hamiltonian Hi in Eq. (H) 

as 

E._ = al + Xal + Xal - - \rl (7) 

where A is a real constant. In this formula, is the interior spatial volume 
energy term. Each interior link carries one unit of energy, and the volume 
energy of a state is proportional to the number of links, splits any interior 
link into two. Since \n)i represents n interior links, af maps \n)i to ?T,|n-|-l) i . 
On the other hand, combines any two adjacent interior links into one. 
Since there are only n — 1 pairs of adjacent interior links in \n)i, a\ maps 
to {n — l)|n)i . Finally, the last two terms represent boundary volume 
energy, but notice that the two boundary links contribute negative energy, 
in total minus one half the energy of an interior link. We stress that this 
value of the relative size of the volume energy contributions is crucial for the 
existence of a continuum limit as will be seen. By now, it should have been 
obvious that ifi, featuring the physics of spatial homogeneity and locality, 
is among the simplest Hamiltonians one can conceive for a spacetime with 
boundaries. 

We now proceed to evaluate the transition amplitudes 
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in the continuum limit. It is convenient to work with its generating function 

oo 



kA=l 



where x, y are complex variables. Introducing 

e{k,i]n) = {i\m\k)i 

and its generating function 

oo 

eix,y;n) = x%^e{k,l; 

k,l=l 

for non-negative integers n, we have 



n) 



Gi{x,y;t) = J2 7-'S>ix,y;n). 



n=0 



(8) 



Note, firstly, from 



Q{k, l-l) = {l\H.\k)^ = [k--]5iu + Xk6i,k+i + Xik-1) 5i,k-i 



that 



eix,y;l) 



xy{l + Xx + Xy) xy 



(1 — xyY 2(1 — xy) 



(9) 



Secondly, since |l)i, |2) i, . . . , and so forth form an orthonormal basis of 71, 
we have 0] 

/dz 1 
-— e(x,-; 1)6(^,1/; (10) 
Zmz z 

where both Q[x^z^^\ 1) and Q{z,y]n — 1) are treated as complex functions 
of z, and the contour encirles all singularities of Q{x, z~^]l) but none of 
Q{z,y;n— 1). Using Eq.(|^), we find that Eq.(|TUp leads to 



Q{x,y;n) 



1 A' 



2 X 



- + - ] + X + x + Xx^] — 



d 



dx 



Q{x,y;n- 1). 
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By Eq.(||), this yields 
d 

—Gi{x,y;t) + 



\2 X J ^ ^ dx 

Together with the initial condition 



Giix,y;t) = 0. (11) 



Gi(x,i/;O) = 0(x,y;O) = -^ (12) 
2 1 — xy 

this first order partial differential equation determines Gi{x,y]0) uniquely. 

We are presently only interested in evaluating the continuum limit of Gi. 
Singularities appear when the coefficients in the square brackets in Eq.(|lTD 
vanish, i.e., for x = y = ±1 and A = =f|. These are identical to the critical 
values found in Refs.|Q and p, and we can apply the same scaling procedure. 
Hence we set 

t = —,x = e-^\ y = e-^\ and A = --e'^', (13) 
a 2 

where a is the distance cutoff and T, X, Y , and A are finite renormalized 

values of t, x, y, and A, respectively, and define the continuum continuum 

limit of Gi, for which we use the same notation, by 

Gi (X, F; T) = lim aGi (x, y; t) . (14) 
Substituting Eq.(|T^ into Eqs.(|TT]) and (|T2]) yields the limiting equation 

^G . (X, Y; T) + (X^ - a) 1 (X, Y- T) + XG i (X, F; T) = 

with the initial condition 

Gi(X,r;0) = — i — . 
2^ ' ' ^ X + Y 

These are identical to the equations found for a model of Lorentzian 
gravity with boundaries in Ref . [p.4[ . By inverse Laplace transformation of 



the solution with respect to X and F, one finds (see Ref.[||) the continuum 
limit of the transition amplitude expressed in terms of the physical lengths 
L = k ■ a, L' = I ■ a of the two boundary components to be 



Gi{L,L'-T) = }/^e-^(^+^')^°*^(^^)/o ( Jy^M^] , (15) 

^ sinh(Tv^) Vsinh(TVA); ^ ^ 

where Jq is the zeroth modified Bessel function. Consequently, our string bit 
model and this particular model of Lorentzian gravity belong to the same 
universality class. We will come back to this model again in Section ^. 
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3 Closed, Non-homogeneous Space-time 



In this section we discuss an example of a Hamiltonian model of cylindrical 
space-time. Conceptually, it is not the simplest such model, which we defer 
to the next section. However, it has the virtue of being solvable by the 
generating function technique of the preceding section, which is our main 
motivation for considering it here. The model is spatially non- homogeneous 
in the sense that the equal time slices have one marked link, that is a link 
created by a matrix of creation operators different from those creating the 
rest. A model of Lorentzian gravity with a marked initial loop has been 
studied in Ref. |Q, and we will find that its continuum limit is reproduced by 
our model. 

The Hamiltonian we consider is given by 

H, = Ti{a^a + b^b) + ^TT{a^a^a+-a^b^b+-b^a^b 

\/ N 2 2 



+ b^ba + b^ab) (16) 

for = oo. As mentioned, this quantum matrix model requires a second 
matrix of creation operators besides a^. The entries of this matrix are written 
as whose corresponding annihilation operator is They satisfy the 
same canonical commutation relations as the a-operators and commute with 
these. 

Let Tc be the Hilbert space spanned by all states of the form 
l^)^ = ]^Tr[6t(atri] 
where n is an arbitrary positive integer. These states form an orthonormal 



basis for in the large- A^ limit p, |T2[ : 



lim (m|ra)c = Smn- 

We consider \n)c as the state of a closed universe with n links, one of which, 
created by b\ is marked, (c.f. the string bit interpretation in Refs.[0 and 

m 

The operators acting on which are relevant to us may be written either 

as 

9f = ]7(^^Tr [(at)^a'] , (17) 



9 



where k and / are positive integers, or as 

1 



]\j{ki+k2+h+l2)/2 



Tr 



(18) 



where ki, /c2, h, and I2 are non-negative integers, (cxf and gf are the restric- 
tions of the same operator to the open and closed string sectors, respectively. 
They are elements of different Lie algebras 0, |13[ , a fact we will see and use 
in Section^.) In the large- limit 



12 



9i\n)c 







and 



fcl,fc2l \ 



\n 



(n-l) \n-l + k)c + 0{ 



ll-l2 + ki + k2)c + OU 



N ' 



if / > n or 
if Z < ra, 



ii li + I2 > n — 1 OT 
if I1 + I2 < n. 



Thus, replaces any adjacent / unmarked links in \n)c with k unmarked 
links, and g^^f^ replaces adjacent I1 + I2 + I links, where the (/i + l)-th link 
is marked, with ki + k2 + l links, where the {ki + l)-th link is marked. Note 
that these operators always preserve the marked link. 

Using Eqs.(|l3) and ([Tsl) , we can paraphrase He in Eq.(^) as 



1 I 0.0 , \ 2 , 1 

9i + ^0,0 + A (^^1 + - 



1,0 , 0,1 

2^^0,0 + 2^0,0 



1 I 0,0 I 0,0 
92 + fl'1,0 + 5-0,1 



where A is a real constant. In this formula, g\ + g^'^ is the volume energy 



term. The terms g\ + l/2gQ'^ + l/2(7o q implement splitting of any unmarked 
link into two unmarked links or splitting the marked link into a marked 
and an unmarked link. Finally, the terms g\ + g^'^ + 5^0,1 combine a pair of 
juxtaposed unmarked links into one unmarked link or combine the marked 
link and a juxtaposed unmarked link into the marked link. Notice that the 
relative constants of the terms in He are chosen such that its action treats 
the marked link in the same way as the unmarked ones. More specifically. 



0,1 



(^l + 5o,o) \n)c = n\n). 



2 , 1,0 I 1 0,1 
9i + 2^0,0 + 2^0,0 



/ 1 I 0,0 I 0.0\ I \ 
[92+ 9i'fl+ 9q:i) \n)c 



= n\n + l)c , 
n\n — l)c. 
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Thus this form of He appears to represent the most natural nearest neigh- 
boring interaction on but notice that it is non-Hermitian. 

As aheady mentioned, it turns out that the model so defined can be solved 
by the same method as that of the preceding section. Since the differences 
between the calculations are only minor, we will skip the details. Using the 
scaling conditions (|l^) one finds that the the continuum limit Gc{X,Y;T), 
defined by the same procedure as in Section 2, fulfills 

Y- T) + (X^ - a) ¥■ T) + 2Xa(X, F ; T) = 

with the initial condition 

a(X,F;0) 



X + Y 



These equations are identical to those found for a Lorentzian gravity model 
in Ref. 1^ with one marked link in the entrance loop. Consequently, the two 
continuum limits coincide, as claimed. For later reference we note that the 
solution in terms of the length variables is |^ 

Ge{L, V- T) = ^;^,-VA(L+L')coth(VAT)^ / Jh/KW_\ 

^ ' ' ^ V i^' sinh(v^r) Vsinh(v^r); ^ ' 

4 Closed, Homogeneous Space-time and Ten- 
sor Product models 

Marking a link in a boundary loop is a convenient technical device in tri- 
angulated models and the relation to the same model with no marking is 
simple, since the marking only gives rise to a factor equal to the length of 
the marked loop in the counting of triangulations. The relation is of a differ- 
ent kind for Hamiltonian models but, as we will immediately see, still quite 
straightforward. 

In spatially homogeneous models, all links in an equal-time slice have 
identical status, meaning that only one type of creation and annihilation 
operators is involved. The simplest nearest neighboring Hamiltonian is then 
given by 



Hi = Tr 



VN 



(20) 
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with N = oo and A a real parameter. Comparing Eqs.(]T^) and (PUj), we see 
that removing the marked hnk restores not only cyclic symmetry but also 
Hermiticity. 

The Hilbert space Ti on which Hi acts is spanned by all states of the 
form 



Tr(at)"|fi), 



where n is a positive integer. In the large- limit P, 12 



N 



lim {m\n)i = n6r, 



Hence |l)i, |2)i, and so on form an orthogonal basis for Ti. We think of 
\n)i as the quantum state of a closed one-dimensional universe topologically 
equivalent to a regular polygon with n sides (c.f. the string bit interpretation 
in Refs. and In terms of the operators introduced in Eq.(|T^, the 

Hamiltonian Hi can be rewritten as 

Hi = gl + Xgl + Xgl 



In the large-A^ limit M, 12 



9l\n)i 

alHi 

^2|l)l 

911^)1 



n\n)i if n > 0, 
n\n + l)i if n > 0, 
0, and 

n\n — l)iif n > 0. 



Here gl is again a volume energy term, gf splits any link into two and g^ 
combines any pair of adjacent links into one. 

It turns out to be tricky, if not impossible, to apply the previously used 
generating function technique to work out the transition amplitude of this 
Hamiltonian in the continuum limit. Instead, we will derive it by diagonal- 
ising H. Before doing so, however, we will digress for a moment and make 
some observations about the underlying Lie algebras of string bit models. 
This will lead to the introduction of a more general and abstract class of 
Hamiltonian models including the one just defined as well as the model in 
Section 2 and, in the continuum limit, the tensor product type models of 
Ref. [0 as special examples. 
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It was shown in Ref. [TS] that if we take Ti as the defining representation 
of the Lie algebra C*i generated by all (7f 's, then under the identifications^ 

gl ^ Lo, g\ ^ Li, and gl ^ 

they satisfy the Lie brackets 

[LqiLi] = —Li, 
[Lo,L„i] = and 

[Li,L_i] = 2Lo, (21) 

which the reader may verify directly and easily. This Lie algebra is nothing 
but the SI2 subalgebra of the Virasoro algebra. Furthermore, since 

gl\l)i = 0, (7}|l)i = |l)iand(l|l)i = l, 

1 1)1 plays the role of the highest weight vector in the defining representation, 
and the highest weight h = 1. It should thus be natural for us to consider 
the model in which the Hamiltonian 

H = Lo + AL_i + ALi 

is an element of SI2 and acts on a certain highest weight represetation. Re- 
call that for general h > the representation space is the Hilbert space 
spanned by the vectors 

|n + l)/. = ^L^i|l),„ (22) 
nl 

where \l)h is the highest weight vector and n is any non-negative integer. 
The actions of L_i, Lq, and Li are given by 

Lo\n + l)h = {n + h)\n + l)hiin>0, 

L_i\n+l)h = {n + l)\n + 2)hifn>0, 

Li\l)h = 0, and 

Li\n + l)h = {n-l + 2h)\n)hiin>0. (23) 



^There is an important difference between the way the Virasoro generators arose in 
Refs.|l^ and and the way they arise here. In those two articles, every Virasoro 
generator L„ was identified as a coset of certain elements of the Lie algebra Ci or Ei, 
a subalgebra of Gi.i. These cosets satisfied the Witt algebra, i.e, the classical Virasoro 
algebra. Therefore, the Witt algebra was a quotient algebra of Ci or Si. Here, on the 
other hand, Lq and Li (and nothing else) are identified with specific elements of Ci 
or a variant of Si. They turn out to form the SI2 subalgebra of the Virasoro algebra. 
Therefore, SI2 is a subalgebra of Ci or the variant of Si. 
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The inner product on is uniquely determined by the commutation relations 
([21|), = 1, Lg = Lq, and L\ = In particular, 

We note that the inner product is positive-definite if and only \ih> ^. 

Next, we revisit the model in Section |^ for a spacetime with boundaries. 
Its Hamiltonian iJi was given by Eq.(|^). The following observation is a 

slight modification of the results concerning the Lie algebra Si generated by 
all (jf's in Ref. [0]: make the identifications^ 



(j\ — -/q ^ Lq, o\ ^ Li, and ^ L_i. 

Then the actions of L^i, Lq, and Li on 71 satisfy the Lie brackets (^). We 
are thus again led to the Hamiltonian H of the form (|l|) . Since 

^2|l)i = 0, (al - l/o) = and = 1, 

|l)i plays the role of a normalised highest weight vector. 

Based on these observations, we will, in the following, consider Hamilto- 
nians of the form (|I]) for arbitrary positive highest weights h. It turns out 
to be possible to diagonalise H for all such h and to evaluate the continuum 
limit of the transition amplitude. Before we perform this task in the next 
section, a few remarks on the interpretation of these models are in order. 

As is well known from the representation theory of s/2, its /i = 1 highest 
weight representation is the symmetric tensor product of two copies of the 
h = 1/2 highest weight representation. This fact has a clear interpretation in 
terms of the gravity models as follows. Taking two copies of the Hamiltonian 
Hi, we define the symmetrised tensor product state 



2 



\ny = \k)i (g)\n - k + l)i 



k=l 



for n > 1 and regard it as representing the state of a closed polygon obtained 
by gluing two polygonal lines of total length n at both ends, where each of the 



ibid. 
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four boundary links of the two polygonal lines has a negative length of —1/4. 
The action of Hi (g) 1 + 1 (g) Hi on the space spanned by the vectors \n) will 
then be easily seen to equal to that of H for h = 1 under the identification 
\n)' = \n). Similar remarks apply to the generators Lq, and Li. 

Corresponding considerations of tensor products and gluing constructions 
for Lorentzian gravity models were discussed in Ref . ^ , as were extensions to 
multiple tensor products. The latter lead to the so-called p-seamed transition 
amplitudes. As will also be seen from the explicit solution in the next section, 
these are reproduced by our algebraic model for integer values of 2h. 



5 Solution to s/2 Gravity Model 

In this section we will show how to obtain the continuum limit of the models 
with Hamiltonian given by Eq.(|T]) for arbitrary h > 0. We will first prove 
that H is diagonalisable and determine the exact energy spectrum. Then 
we will determine the asymptotic form of the eigenvectors close to criticality, 
which will enable us to extract the continuum limit. 

As already remarked, Th is spanned by the vectors \l)h, |2)/i, and 
so forth defined by Eq.(^). These vectors being orthogonal, it follows that 
states in are given by 

oo 

^an\n)h, (25) 

n=l 

where 

oo 

J2n'"'-^K\'' <oo (26) 

n=l 

by Eq. (P^) . Clearly, H is an unbounded operator defined, e.g., on vectors for 
which the sequence is rapidly decreasing. 



5.1 Diagonalisation of H 

We will apply a refined version of the Frobenius method used to solve a very 
similar quantum matrix model in Ref.]l2[. Assume, for E > 0, that 

oo 

(t>=Yl (^n\n)h (27) 

n=l 
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is an eigenstate of H. Using Eqs.(0) and (P^, we may write the eigenvalue 
equation 

H<p = E<p 

as 

\{n + 2/i)a„_|_2 + {n + h — E)an+i + Xnan = 

for all non-negative values of n. (The value of the new unknown oq is imma- 
terial because its coefficient is 0.) Asymptotically for large n, the equation 
reduces to 

Xan+2 + an+i + Aa„ ~ 0, 
whose solutions are of the form 

an ^ ap" + 

where 

in the continuum limit. Since \p\ < 1 for |A| < 1/2, it follows that a„ must 
asymptotically behave as p"^ in order that be normalisable in accordance 
with Eq. (|26|) . Hence, we set 

an = hnp"", 

resulting in the equation 

A(n + 2h)hn+2P^ + {n + h — E)bn+iP + Xnbn = 

for all non-negative values of n. 
In terms of the increments 

A6„ = bn+i - bn and A^6„ = A6„,+i - A6„, 

this equation may be rewritten as 

\{n + 2h)pA%n + [(2Aj9 + l)n + {AhXp + h - E)] Abn 
+ {2hXp + h - E)bn = 0. 



Introduce the ansatz 15 



K = Y1 

r=0 

16 



(r) 



where the real numbers Cr depend not on n but on r only, and the factorial 
polynomial n^^^ is defined by 



(r) 



n{n — 1) • ■ ■ (n — r + 1) if r > 0; and 
1 if r = . 



From 

oo oo 

Abn = ^ rCrn^^~^^ and A^6„ = ^ r(r — l)c, 

r=0 r=Q 

we then obtain the equation 

\p{r + 2)(r + l)(r + 2h)cr+2 + [(3r + Ah)\p +{r + h-E)]{r + l)c^+i 

+ [(1 + 2Ap) {r + h)-E]cr = (29) 

as a sufficient condition for the coefficients Cr to fulfill for all values of r. 
For a non-negative integer R, we now set 

E = Er= (1 + 2Xp) {R + h) (30) 

and see that we obtain a unique solution for c^, up to constant multiples, for 
which = for r > R. The corresponding eigenstate is thus of the form 

oo 

= Ilc'i?('^KI'^)/., 

n=l 

where Cr^ti) is a polynomial of degree R in n. 

Since H is Hermitian, the found eigenstates form an orthogonal set. More- 
over, it is a complete set, which can be seen as follows. Take any vector in 
Tfi given by (p5|) and assume it is orthogonal to all eigenvectors Since 
the Cr(?t,)'s span all polynomials, this means that 

n\r{2h) 

where S is an arbitrary non-negative integer and we have used Eq. (|2^) . Mul- 
tiplying this equation by / S\ and summing over S give 

^ ^T{n + 2h) 
""^ n\r{2h) 
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for \z\ < — \ogp. Obviously, the left hand side is an analytic function of z 
in the half plane '^{z) < — \ogp and hence vanishes there. Restricting z to 
the imaginary axis, we obtain a vanishing Fourier series, and consequently 
its Fourier coefficients vanish. This proves that = for all n > 1 and the 
completeness of (pn for all i? > follows. 

Thus Eg. (pOD gives the whole energy spectrum. Note that in the limit 

(H), 

as a ^ 0, so the model is well behaved in this limit. 



5.2 Asymptotic behaviour of eigenstates 

In order to determine the continuum limit of the transition amplitude we will 
need the asymptotic behaviour of the polynomials Cii{n) under the scaling 
conditions given in Eg. ([T3|) . Since n scales as we need to exhibit the 
leading behaviour of the coefficients in C^ln) as a — > 0. Make the ansatz 
that Cr+i and aCr are of the same order in the small-a limit. The recursion 
relation (p9|) then gives 

R-r 



Cr+l 



2\fXaCr. 



(r + l)(r + 2/i) 
Iterating this eguation yields 

-2^aYV{2h) ( R 



co- 



r(r + 2h) 



r 



which, owing to the scaling oi n = L ■ a ^, yields the asymptotic form 

where all summands are of order 1. 

The behavior of cq is fixed by reguiring that ipR be normalised. Using the 
fact that ipR is orthogonal to all vectors in Th of the form ( p5l) with = n'^p", 
for s = 0, 1, . . . , and i? — 1, we get 

oo oo ij 
n=l m=l r=l 

- |C«(")CH«V"M. (32) 
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Substituting Eq.(PTD into Eq.(^) results in 



-2 



,t^or(i? + 2/i)r(r + 2/.) V ^ n! 
^ (-2VAa)^+T(2/i) / i? \ r(n + i? + r + 2/i) 3^ 

where, in the last step, we have used that for k real 

r(n + k) 



1 as n — >■ cxo 



By the binomial theorem and the relation 

~ 1 - 2VAa, 

we finally obtain 



c, 



-2 



:^ (-2v^a)^+"r(2/i) / i? \ r(i? + r + 2/1) 



° Sr(^ + 2/i)r(r + 2/i) \^ r y (l_p2)R+r+2h 

^ 1 Vf ^ »+, r(i? + r + 2/^)r(2/i) 

(2VAa)2'^i^o r(i? + 2/i)r(r + 2/i) 

1 RW{2K) 
{2VAay^r{R + 2h) ' 

where, in the last step, we have made use of the identity 

f^^^ ^ WT{r + 2h) \ r ) 
which is a special case of the Chu-Vandermonde identity. (See, e.g., Ref.[]I^].) 

5.3 The continuum limit 

We are now ready to compute the continuum limit of the transition ampli- 
tude. The unnormalized transition amplitude is defined by 

GuiL,L';T) = lima"(-|e-*^|-),, (34) 

a-*o a a 



R 

r 



(33) 
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where t and A are given as in Eq.(|TB|), and the exponent a is to be determined 
such that the hmit exists. Recall from Eq.(p4D that the states \^)h are not 
normalised. On the other hand, the more natural amplitude defined in terms 
of the normalised states may simply be obtained from Gu by Eg. p^) ; we will 
come back to it later. 

Inserting the complete set of states {4>r}, we have 

OO Tl T 

G^{L,L';T) = \ima''Y.^~'''H-\<l>R)h{(l>R\-)h 

OO T I T 

^ lima" E e-2^(^+^)(- | | . (35) 

Using Eqs. (PT|) and ( P^D as well as 



we have 



JR —/h 

a 



\a / a a 
(2Va) 



T{R + 2h) ^2h-i j2 f R\ir^ 



\ mT{2h) f^^ r(r + 2h)\r ^ 

Inserting this expression into Eq.(p5D and choosing 

a = 2/1-2, 

we find that the continuum transition amplitude exists and takes the form 
GuiL, L'; T) = 

^ > ^ > Rir{2h) ^,tonr + 2hns + 2h) 

' ( ^ ) ( ^ ) ^"^^■^^■^'^e"^(''+^)^^. (36) 

A priori, it is not obvious that this series is convergent for all positive 
values of L ,L', and T. One way to see this, and simultaneously obtaining a 
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more manageable expression for G^) is to apply an integral representation of 
the reciprocal Gamma function (see, e.g., Ref.||T7[|): 



Tix) 



[z + hy 



(37) 



In this formula, x and j3 are arbitrary positive numbers; e is real and h 
complex, and they satisfy 3?(6) > e; and the branch cut of {z + 6)^ lies on the 
negative real axis if x is not an integer. Apply Eq.(|57D to l/r(r + 2h) with 

x = r + 2h, (3 = VAL, b = l,z = X, 

and e positive and infinitesimally small, and to l/r(s + 2h) with 

x = s + 2h, (3 = VAL', 6 = 1, z = Y, 

and the same e in Eq. ( |5B| ) . The binomial theorem can then be used to perform 
the summation over r and s. Apply once more the binomial theorem to the 
sum over R, and keep the branch cut of every [■ ■ ■]^'* on the negative real 
axis. With this choice of the branch cuts, {ziZ2Y^ = z\^Z2^ if both '^{zi) and 
'Si{z2) are positive. Using this fact, we finally get 



Gu{L,L';T) = 4'^A 



h A 1—h 



27Ti 



€+ioo 



27ri Je- 



[Y + 1-{Y- l)e~^VATj2h 



X 



y+i+(y-i)e-2^^ 

Y+l~{Y~l)e-'2^T 



2h ' 



Integrals of this type have been carried out for integer values of 2/i in 0, 
but can be obtained for arbitrary values of 2h. Indeed, apply Eq. ( ^71) to the 
integration with respect to X and we immediately get 

Ah„-2h^T.nrT2h~l 



e+ioo 6Xp 

dV 



T{2h)2m 



y(i-e 



-2-/AT 
-2-/AT 



[¥{1 



,-2 



v^^) + (1 + e-2v^^)] 



2h 



Inserting the expansion 
exp 



/Tr y(l+e-^^^)+(l-e~^v^T) 



exp 



Val 



l+e- 



-2'/AT 



ra=0 



l_e-2v^T Y{l-e- 



2 VAT 



_L 

)+{l+e 



-2\/AT 



l_e-2%/AT 
n 
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applying Eg. P7|) again to the remaining integration, and recalling that 

oo ^2n+2/i-l 

hh-i{z) = 22 i ^• 

^on!r(2/i + n) 

is the {2h — l)-th modified Bessel function [113, we finally obtain 

^ ^ r(2/i) sinh( VAT) \smh{VAT)J 

Define the normalized transition function G{L, L'; T) by normalizing \ -)h 
in Eg. (p4D . It follows from Eg.(p^) that we have to choose a = 1 and that 
G{L,L';T) deviates from Gu{L,L';T) by a factor of {LL'f-^ /T{2h). Con- 
seguently, 



We note that if we put h = 1/2 in Eg.(^), we will obtain Eg.(|15]); for 



h = 1 we obtain Eg.(^) except for the factor y L/L', which originates from 
the non-Hermiticity of He, if h is an integer, Eg.(^) will coincide with the 
propagator calculated in the proper-time gauge of two-dimensional guantum 
gravity in Ref.||TT[| with the winding number h — 1; finally, if 2h is an integer, 
Eg.(R3) will be exactly the (2/i)-seamed correlation function in Ref.[H. 



6 Discussion 

We have, in this paper, investigated different Hamiltonian models of two- 
dimensional guantum gravity. Clearly, one open problem is the proper phys- 
ical interpretation for the s/2 gravity model with a non-integer value of 2h. 
Perhaps this describe the interaction of gravity with matter, an important fu- 
ture problem on its own right. One could, for instance, study one-dimensional 
guantum spin systems whose Hamiltonians couple spin configurations of dif- 
ferent sizes. Diagonalisation of such Hamiltonians seems to pose interesting 
new problems. Another possibility is to extend the class of Hamiltonians 
defined in this paper by exploiting the representation theory of the full Vi- 
rasoro algebra instead of the s/2 subalgebra; such models might describe the 
coupling between matter and gravity. 
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Finally, extension to higher dimensional Hamiltonian models of quantum 
gravity is an ultimate goal. It is rather straightforward to produce candidate 
models of nearest neighbouring type, but extracting non-trivial information 



from such models seems a non-trivial task. We refer to Ref. [0 for a recent 
work on higher dimensional Lorentzian gravity. 
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